Spin glass behavior of frustrated 2-D Penrose lattice in the classical planar model 
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Via extensive Monte Carlo studies we show that the frustrated XY Hamiltonian on a 2-D Penrose lattice 
admits of a spin glass phase at low temperature. Studies of the Edwards- Anderson order parameter, spin 
glass susceptibility, and local (linear) susceptibility point unequivocally to a paramagnetic to spin glass 
transition as the temperature is lowered. Specific heat shows a rounded peak at a temperature above the 
spin glass transition temperature, as is commonly observed in spin glasses. Our results strongly suggest that 
the critical point exponents are the same as obtained by Bhatt and Young in the ± J XY model on a square 
lattice. However, unlike in the latter case, the critical temperature is clearly finite (nonzero). The results 
imply that a quasiperiodic 2-D array of superconducting grains in a suitably chosen transverse magnetic 
field should behave as a superconducting glass at low temperature. 
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The notion of frustration, first introduced by 
Toulouse 1 , has played a key role in the theoretical un- 
derstanding of spin glasses. It is now believed that frus- 
tration is a necessary condition for the existence of a 
spin glass phase. However, all real systems (alloys) ex- 
hibiting spin glass behavior have an additional common 
feature, namely the topological disorder of magnetic com- 
ponent^). On the theoretical front, models of uniformly 
frustrated spin systems without disorder have often failed 
to show spin glass behavior, and the debate as to whether 
disorder, in addition to frustration, is necessary in a spin 
glass continues. In this paper we present strong evidence 
that a frustrated XY Hamiltonian on a 2-D quasiperiodic 
lattice admits of a low temperature spin glass phase. The 
degree of frustration in this case varies widely from one 
lattice point to another, but in a predictable way deter- 
mined by the quasiperiodic structure. 

We consider the Hamiltonian for the XY model de- 
scribing the interaction between 2-D spin vectors with 
orientations 9i and 0j situated at lattice sites i and j via 
a nearest neighbor coupling parameter J: 



- 0i + A t . 



(1) 



where the summation is restricted to nearest neighbor 
pairs [ij]. The parameter Ay controls the frustration in 
the model. In the context of an array of superconducting 
grains, Qi is the phase of the grain i and the above Hamil- 
tonian can be seen as describing the resulting Josephson 
junction of the grains "minimally coupled" to a trans- 
verse magnetic field with vector potential A with 
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$o is the elementary flux quantum ^ associated with 
the Cooper pairs, and r\ denotes the lattice sites. Here 
the magnetic field acts as the source of frustration: an 
Aij which is an odd multiple of n essentially renders the 
bond [ij] negative. 



The directed sum of A^ about a plaquette in a 2- 
D lattice can be written as 2wf, where / is the flux 
through the plaquette in units of <I>o. The 2-D Pen- 
rose lattice 2 is composed of two (fat and thin) rhombic 
unit cells (plaquettes). Since the ratio of the areas of the 
fat and the thin rhombuses in the Penrose lattice is the 
Golden Mean (r) , only one set of plaquettes can be fully 
frustrated at a time with a suitable choice of the magnetic 
field giving / = 1/2. The flux / through the individual 
plaquettes in the other set will then be an irrational num- 
ber. 

In what follows we will present results for the case 
where the 'thin rhombus' is fully frustrated, mentioning 
at the outset that the results for the fully frustrated 'fat 
rhombus' case are qualitatively similar. All our results 
are obtained via Monte Carlo (MC) simulation based on 
the Metropolis algorithm 3 using periodic boundary con- 
ditions on periodic Penrose lattices. We have cooled our 
systems in a quasi-static manner, starting from a high 
temperature (T(in units of J)> 2.0) random configura- 
tion and then heated the system in the same quasi-static 
fashion. Since we performed the simulation in n blocks, 
the heating and cooling data are obtained by averaging 
over these blocks, with the error bars representing the 
standard deviation, obtained by dividing the square root 
of sum of squares of the deviations from the mean by 
y/n—1, instead of y/n. We then perform a 'grand aver- 
age' over the heating and cooling data. 

Magnetic moment per lattice point calculated for peri- 
odic Penrose lattices of various sizes is found to be small 
(< 0.02) over the entire temperature range. The mag- 
nitude of the moment decreases steadily with the size of 
the cluster, suggesting that the magnetization is strictly 
zero in the thermodynamic limit. To study the freezing 
of the spins at the lattice sites we calculate the Edwards- 
Anderson 4 order parameter, a measure of breakdown of 
ergodicity in the system. In Fig.l we show this order 
parameter defined by 



Qea 



(3) 
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where Si — (cos $i, sin 9i), and ( ) denotes a canonical 
ensemble average. In a completely frozen system qea is 
unity while for a completely ergodic system it is zero. 



N = 4414 
L, = 64.9 
L y = 55.2 



0.4 o.e 
kJ/J 



FIG. 1. Edwards- Anderson order parameter as a function 
of temperature for N=4414. The linear dimensions of the 
lattice in x and y directions are indicated as Lx and L y , 
respectively 

This order parameter shows a monotonic decrease with 
increasing temperature, clearly vanishing at tempera- 
tures beyond 0.5 for the 4,414 site periodic Penrose lat- 
tice shown in Fig.l. These results were obtained by av- 
eraging over 5 blocks of 60,000 configurations, generated 
after equilibrium was achieved. Since the vanishing of 
the order parameter with a long tail is a consequence 
of the finite system size, it is expected that the tail re- 
gion will decrease with increasing system size and even- 
tually disappear in the thermodynamic limit. However, 
we find that this tail persists, even for our largest system 
(N = 11,556) and, consequently, the spin glass tran- 
sition temperature cannot be appropriately determined 
from Fig.l. Thus, we use other quantities to provide an 
estimate of the transition temperature Tf. It is clear, 
however, that the system has a low temperature phase 
with a nonzero order parameter qEA- 

Local (linear) susceptibility per spin, calculated from 
the fluctuations in the magnetization (net magnetic mo- 
ment m for a lattice of N sites), 
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125,000 MC steps, while 5 blocks of 15,000-45,000 steps 
were used for higher temperatures. The large hystere- 
sis in the low temperature region indicates a high num- 
ber of metastable states, which is a characteristic of spin 
glasses. These metastable states give rise to large error 
bars in (m 2 ) — (to) 2 at low temperatures, which are fur- 
ther accentuated by a division by T in Eq.(4). Although 
we feel that it might be possible to reduce the size of 
these error bars, this would require very long runs and 
one must also ensure that the system does not become 
trapped in one of these metastable states. Nevertheless, 
despite the large error bars, a cusp- like feature in x is 
clearly visible at Tf ~ 0.15. We find a saturation in 
this cusp with respect to system size, which is consistent 
with spin glass behavior. Note that the high tempera- 
ture (above Tf) phase is strictly paramagnetic. Here the 
decrease in \ is Curie-like, with T\ becoming a constant 
as shown in the inset of Fig. 2. 
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is shown in Fig. 2. At low temperatures, (0.02 — 0.2) 
the results are obtained by averaging over 5 blocks of 
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FIG. 2. Linear susceptibility as a function of temperature, 
with a peak around T=0.15. 

Specific heat obtained from the fluctuations in the en- 
ergy is shown in Fig. 3 for various system sizes. The re- 
sults are averages between heating and cooling, with the 
low temperature results having a somewhat larger hys- 
teresis. All averages were obtained after equilibriating, 
however, the high temperature values were obtained by 
averaging over 5 blocks of 15,000 steps, whereas 5 blocks 
of 45,000 steps were used for the low temperatures. The 
result for the unfrustrated case is shown in the inset. 
Both frustrated and unfrustrated cases show saturation 
in specific heat with respect to system size. In the unfrus- 
trated case the peak in the specific heat is at 1.10, which 
is beyond the Kosterlitz-Thouless (KT) transition 5 tem- 
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perature Tkt- The specific heat peak for the frustrated 
case is more rounded relative to the unfrustrated case 
and occurs at a temperature lower than Tkt- This tem- 
perature is, however, higher than the temperature Tf at 
which, we believe, the Edwards-Anderson order parame- 
ter goes to zero or a cusp appears in the linear suscepti- 
bility. 
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FIG. 3. Specific heat as a function of temperature. The 
inset shows the specific heat for the unfrustrated case with 
the peak slightly beyond the KT transition temperature 

The saturation in the peak height of the specific heat is 
a consequence of the fact that it appears at a temperature 
at which the spin glass correlation is finite. Note that in 
both the unfrustrated and the frustrated cases the zero 
temperature specific heat approaches a value of 0.5ke, 
consistent with a linear spin- wave theory. 

In a ferromagnet, the approach to the ferromagnetic 
phase from temperatures above the Curie temperature 
Tc is accompanied by a dramatic increase in the range 
of the spin correlations, which then diverges at Tc- A 
corresponding phenomenon occurs in spin glasses. How- 
ever, it is not the spin correlation function (Si ■ Sj), but 
rather its square that acquires a long-range. This leads to 
the divergence, at the spin glass transition temperature 
Tf, of the spin glass susceptibility 



(T > T f ) 



(5) 



X SG satisfies a finite-size scaling relation of the form 5 

XsG =L^x{L lJv {T-T f )), (6) 



where x is the scaling function, L is the system length, 
v is the exponent for spin glass correlation length £ for 
T > Tf, and 77 describes the power law decay of the spin 
glass correlation at Tf. For a square lattice L = y/~N, and 
the above scaling relation can also be expressed in terms 
of N. We have used this scaling relation, expressed in 
terms of N, for periodic Penrose lattices of various size 
to study Xsg ■ To ensure a proper convergence of x S g > 
we have averaged over 5 blocks of 40,000-60,000 steps at 
low temperatures (< 0.2) and and 5 blocks of 60,000- 
80,000 steps at higher temperatures. By examining our 
results every 5,000 steps, we find little change in x S c 
over the last 5,000-10,000 steps. Thus, we estimate that 
these chain lengths produced at least a 95% convergence 
in Xsg- 
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FIG. 4. Scaling plot for spin glass susceptibility. The inset 
shows the log-log plot of the Edwards- Anderson order param- 
eter with (T-T/) with T/ =0.137 obtained from the fit to the 
scaling relation (6) with 77=0.02 and u—2.6 , as given by Bhatt 
and Young 5 . The slope of this plot yields a value of 0.26 for 
the exponent /3, in agreement with the hyperscaling relation 
(5 = vti/2. 

In Fig. 4, we show the scaling behavior of Xsg m terms 
of N for Tf = 0.137, where we have used the values of 
v = 2.6 and rj = 0.2 reported by Bhatt and Young 6 in 
their study of ± J XY model on square lattices. These re- 
sults represent the average over cooling and heating and 
lie within the 0.02 — 0.45 temperature range. A choice 
of Tf — results in four distinctly separate curves, in- 
dicating a clear breakdown of the scaling relation. If we 
try to obtain a good fit to the scaling curve with Tf = 0, 
then the exponents v and r\ deviate drastically from the 
values obtained by Bhatt and Young 6 and others (see ref- 
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erences in Bhatt and Young), obtained for 2-D XY spin 
glass models. 

Furthermore, we have examined the Edwards- 
Anderson parameter in the vicinity of Tf, where the re- 
lation qEA ~ (T — Tf)P is supposed to hold. We find that 
our value of T/ = 0.137 yields a value of (3 = 0.26, (see 
inset of Fig. 4) which satisfies the hyperscaling relation 
(3 = ^rp valid for our 2-D system. Note that this value of 
Tf is clearly consistent with our results for the suscepti- 
bility which shows a cusp at T ~ 0.15 (Fig. 2). The 
inescapable conclusion seems to be that the frustrated 
2-D Penrose lattice has a spin glass transition tempera- 
ture above zero with exponents that are equal (similar) to 
those obtained for random XY models on square lattices. 

In summary, we have shown, via the Edwards- 
Anderson order parameter, spin glass susceptibility and 
local (linear) susceptibility, the existence of a low temper- 
ature spin glass phase in frustrated 2-D Penrose lattice. 
Our results for magnetization and specific heat also sup- 
port this picture. Our detailed study of the unfrustratcd 
case (not reported here) reveal a behavior similar to that 
of a square lattice 7 with a somewhat higher KT tran- 
sition temperature. We would like to add that similar 
studies by us performed on other frustrated quasiperi- 
odic lattices (such as the octagonal lattice) also show a 
low temperature spin glass phase. 

A few comments are in order at this stage. Our results 
are consistent with those of Halsey 8 , who finds a spin 
glass phase for the frustrated XY model on a square lat- 
tice with an irrational flux through the plaquettcs. Note 
that we can fully frustrate only one of the two elementary 
plaquettes at one time, the corresponding flux through 
the other plaquette being irrational. In both cases, the 
model studied by Halsey 8 and the frustrated Penrose lat- 
tice, the parameter cos(-Ajj), which determines the na- 
ture of the coupling along the bond [ij], shows a similar 
distribution in its values, varying widely between 1 and 
-1. This is different from most of the fully frustrated reg- 
ular lattices studied so far. For example, the fully frus- 
trated square lattice studied by Teitel and Jayaprakash 9 
or the triangular lattice studied by Shih and Stroud 10 , 
where the nature of the transition is seen to change from 
KT to the Ising type, has a much narrower distribution 
of cos( Aij). 

The experimental implication of our study is that a 
quasiperiodic array (such as the Penrose or octagonal 
lattice) of Josephson junctions in a suitably chosen trans- 
verse magnetic field should behave as a superconducting 
glass at low temperature. Via advanced microfabrica- 
tion techniques 11 , it should be possible to generate such 
quasiperiodic arrays of superconducting grains. Experi- 
mental investigation on a 2-D fractal (Sierpinski-gasket) 
network has been reported 12 . Halsey 8 has pointed out 
that for superconducting arrays with low normal-state 
resistivities the glass transition should basically appear 
as a mean-field transition, with fluctuation effects be- 
ing barely observable. For arrays with high normal- 
state resistivities the fluctuation effects will cause the 



glass transition to deviate substantially from a mean- 
field transition, with noticeable system-dependent de- 
tails. An important feature of such glassy superconduc- 
tors, as pointed out by Ebner and Stroud 13 , is a large 
difference between their dc and ac susceptibilities. 

Finally, a comment about the lower critical dimension 
for the spin glass phase. It is widely believed that the 
lower critical dimension for the spin glass phase is greater 
than two 14 . This is corroborated by Monte Carlo studies 
by Binder and co-workers 15 (see also ref. 14) and Young 
and co-workers 6 on square lattices. Our results indicate 
that quasiperiodic 2-D lattices may provide an exception. 
The glassy behaviour we are reporting is truly the behav- 
ior of macroscopic systems and not a reflection of some 
other model masked by the finite system size. This be- 
havior is guaranteed for the particular choice of frustra- 
tion considered in this paper, dictated by the quasiperi- 
odic lattice structure and the magnetic field. Unlike in 
the above studies 6,14,15 with completely random distri- 
butions of the coupling constants along the bonds, the 
effective coupling constants in our model are nonrandom. 
However, our results obtained for the particular magnetic 
field considered are amenable to experimental verification 
and merit attention. 

Financial support for this work was provided by the 
Natural Sciences and Engineering Research Council of 
Canada. 
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The authors regret that the following errors remained un- 
detected during the proof stage. 

(1) In the abstract ! ±J XY model' should read '± J Ising 
model'. 

(2) Rcf. 6 should read: 

(a)R.N. Bhatt and A.P. Young, Phys. Rev. B 37, 5606 
(1988); (b) S. Jain and A.P. Young, J. Phys. C: Solid 
State Phys. 19, 3913 (1986). 

(3) In the caption of Fig. 4, Ref. 5 should be changed to 
Ref. 6(a), and r\ = 0.02 should be changed to r\ = 0.2. 

(4) Ref.5 immediately preceding Eq.(6) should be 
changed to Ref. 6(a). 

(5) Regarding the exponents, note that our results agree 
much better with those of Bhatt and Young 6 ' ) (± J Ising 
model) than with those of Jain and Young 6 ' 6 ) (± J XY 
model). Accordingly, references to ±J XY model study 
by Bhatt and Young 6 should be changed to ±J Ising 
model study by Bhatt and Young 6 '") at appropriate 
places in the main body of the paper. 

One final comment: It is not clear to us at this stage 
whether the agreement with the ± J Ising model is purely 
coincidental, or there is a connection with the result of 
S. Teitel and C. Jayaprakash, Phys. Rev. B 27, 598 
(1983) (Ref. 9) on square lattice, where the nature of 
the transition is seen to change from KT to the Ising 
type as a result of frustration. 

The authors are thankful to A.P. Young for pointing 
out the referencing error regarding the ± J XY and ± J 
Ising model calculations. 
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